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MTHSGASIN ; G Floating Point Si i - mq 
vai pa Point om rovtig MEE FLSA BANMTE DSIRE ess ch] 


-TITLE MTHSGASIN 3 G Floating Point Sine routine 
: (GASIN, GASIND) 
. IDENT /1-005/ : File: MTHGASIN.MAR EDIT: RNH1005 


SERRA EAA A AAA AAAE AAA AAAAAAAAAAAAAAAAAAAAEAAAAREEAAEAAAEERAEEAAAEEEEETEEAAEEE 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
heey eo iy eGORPORATION, MAYNARD, MASSACHUSETTS. 


THIS SOFTWARE IS ora UNDER A LICENSE AND MAY BE USED AND gd LA 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE ay ney oe NOTICE. THIS SOFT ARE OR ANY OTHER 
COPIES hay al oy NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 


DOQOQOQCOOOGOOOOCOO 
ooooo 


ht 8 3 ss 


® 
® 
7 
® 
® 
* 
® 
® 
® 
t 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
® 
® 
® 
te 
* 
® 
* 
® 
2 
& 
* 


. . . . . 
eeeeeeeeeeeeeeeenennenns 


4 : 
+ 
0 : ; 
0 3 
i: 
0 : 5 
44 72 
sit g : TRANSFERRED. 
0900 3 : THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
000 ; ND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
944 y H CORPORATION. 
itd ¢ : DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
4 7: SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
0900 i : RERAATAAATATAAAAEKSKERAAAAAAAAAEAAKARERAARAAAEAAATAAARAEAAAAAARAARAAERAAERAARRAREAREEARE 
0000 8; 
0000 7 3 FACILITY: MATH LIBRARY 
0000 0 ;+ 
44 1 3 ABSTRACT: 
0000 f ! MTHSGASIN is a function which returns the G Hheor ing point arcsine in 
0000 4 ; radians of its G floating point argument. The call is standard 
it 2 : call-by-reference. 
0000 7 : MTHSGASIND is a function which returns the G fleet tne point arcsine 
0000 3 : in degrees of its G floating point argument. The call is standard 
0000 3 call-by-reference. 
000 40 ; 
00 $5 : 
00 rf ; VERSION: 1 
50 45 : HISTORY: 
0 rt 3 AUTHOR: : 
? : Steven B. Lionel, 15-Jan-79: Version 1 
0 i : MODIFIED BY: 
0 0; 
00 13 
00 2; 


. ¢ 6 Ape Bde % MT 
ut Stony’! Deeadied eurrent Fare History ‘'Sesepalope $1:94:49 PAYANTS SReSRrHeRcOG maps1 Page 3, 2 
: -SBTTL HISTORY ; Detailed Current Edit History 
0 
0 
0 


; ALGORITHMIC DIFFERENCES FROM FP=11/C ROUTINE: none 


¢ - Change JSB entry to MTHSGASIN RBG 
= Added degree entry points. RAH 29-MAR- 

4 = Changed shared external references to G* RNH 25-Sept-1981 
005 = Eliminate symbolic short Literals. RNH 15-Oct-81 


-Sept-1979 
gi 
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AOAAAASO 


1- 4 - Adapted from MTHSDASIN yersten {-001, L 15-Jan-79 
1- 6 B 

1-00 

1 

1 


0 
0 
000 
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MTHSGASIN ; G Floating Point Sine routine =SEP-1984 01:24:49 VAX/VMS Macro V04-00 Page aT 
1-008 BECLARATIONS 3; Declarative Part of Modul panies 7 9138 :19 CMTHRTL.SRCIMTHGASIN.MAR; 1 . (3) e~( 
67 »SBTTL DECLARATIONS 3; Declarative Part of Module 
; INCLUDE FILES: OTSPARAMS..MAR 
ie 3 EXTERNAL SYMBOLS: 
te -DOSABL GBL 3; Force error for undefineds 
0 75 -EXTRN MTHSGSQRT_R ; Square root routine 
0 6 eEXTRN MTHSGATAN_R 3; Arctangent routine 
00 4 -EXTRN MTHSGATAND_R7 3 Arctangent routine 
38 4 -EXTRN MTHSSSIGNAC 3 Math signal routine 
+ 4 eEXTRN MTHSK_INVARGMAT ; Error code 
44 13 
000 ¢ 3; EQUATED SYMBOLS: 
448 
00000004 0000 4 value = 4 ; value.rg.r 
0000 85 
0000 a6 3 
0000 7 ; MACROS: none 
0000 H+: 3 
800 3 ; PSECT DECLARATIONS: 
00000000 91 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
44 38 3; program section for math routines 
0000 94 ; OWN STORAGE: none 
0000 95 ; 
0000 96 ; 
8008 97 ; CONSTANTS: 
000 98 ; 
0000 99 
4 100 G_PI_OVER_2: 
2018 5444 21FB 4019 000 101 -@ORD %0040031, *0020773, *0052104, *0026430 
0008 Ie 3; PI/2 
43 103 G_PI: 
2018 5444 21FB 4029 is ise -WORD *0040051, *0020773, Neonat ft *0026430 
0010 198 6_90: 
00000000 80004076 oa t3 4 -LONG “%X80004076, “x0 : 90 


————+ 
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GF eeite, casin  'SoSEp=1986 $1:53;19 EMTMNTESREIRTNGASGN.MARST 2% (a) 2=( 


-SBTTL MTHSGASIN = Standard G Floating GASIN 


MTHSGASIN : 
1-008 , 


FUNCTIONAL DESCRIPTION: 
; GASIN - G floating point function 
; ; GASIN(X) is computed as: 


ATAN(X/SQRT(1=X*##2)), 


; ; CALLING SEQUENCE: 
gasin.wg.v = MTHSGASIN(x.rg.r) 
S INPUT PARAMETERS: 


LONG = 4 3; define eapeary multiplier 
x = 1 * LONG 3; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: G floating arcsine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signals: ae INVARGMAT if {X! > 1 with rygor ved operand in RO/R1 
(copied to the signal mechanism vector CHFSL_MCH_RO/R1 by CISSSIGNAL). 


Associated message is: ‘INVALID ARGUMENT". banal g ~ reserved 
operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 


: IN(X 
: If X = 1. then GASIN(X 
: SINC 
; nG 


00000004 
00000004 


NOTE: This procedure disables floating point underflow, enables integer 
over filow. 


-ENTRY MTHSGASIN, “M<IV, R2, R3, R4, RS, Rb, R7> 

3 standard cal l-byoreterence en entry 

3; disable DV (and FU) e 
MTHSFLAG_JACKET 3 flag that this is a jacket 1 Ben in 


MOVAB G*MTHSSJACKET_HND, (FP) 


PDEA BB BB BE EB EAA AI IWIPPUNINININPINININ 2 2 OO Oo oe “7 


=O OONA NEW 9 OD NAME WIN 9 OD NIA UE WIN 9 OD NAME WINS OO ONOAOU EO 


mee a kk a ed ot od = — — — 9 9 = 9 = = = = SS SS = 3 SS SY YH SY I YS SS SSS 


6D 00000000°GF 9 


set handler address to jacket 
handler 


$$$ $$$} rc 
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MTHSGASIN ; G Floating Point Sine routine 16-SEP-1984 01:24:49 YAX/VMS Macro V04-00 Pp 5 MT 
atts WTASGASIN 2 Standara'e Floating GASIN -SEp-19B6 91:53:19 EMTURTE. SaedaTHGASiN.maR:? 28° (2) 2= 
7° § 3; case of an error in routine 
} ! t ; at on error convert signal to user PC 
: and resigna 
50 04 BC SOFD 0021 165 MOVG @value(AP), RO ; RO/R1 = {Xi = @value(AP) 
01 10 0026 166 BSBB MTHSGASIN_R7 ; call special GASIN routine 
04 Bose 19 RET 3 return with result in RO/R1 


G 6 
MTHSGASIN ; G Floating Point Sine routine 16-SEP-1984 01:24:49 VAX/VMS Macro v04-00 Pa 
vt MTHS RY = Special GASIN routine pat ioet yy 94:33:48 MTHRTL.SRCIMTHGASIN.MAR; 1 - 


-SBTTL MTHSGASIN_R7 = Special GASIN routine 
; Special GASIN = used by the standard routine and direct JSB call. 


; 0< ixi < 1.0 


1 
1 
1 ; 
174 ; CALLING SEQUENCE: 
175 ; save anything needed in RO:R7 
176 ; V TTT a. : input in RO/R1 
177 ; JSB MTHSGASIN_R7 
: : RSB 3 return with result in RO/R1 
1 oO 
181 MTHSGASIN_R?:: 3; special GASIN routine 
182 MTHSGASIN R9:: 3; Release 1 name 
56 39 50F 1 MOVG RO, R6 3 save X in R6/R7 
4 1 BEQL RETURN : return DASIN(O) = 0 
186 ; 
1 3 0< ix! 
: 3 
50 8000 HA AA 1 BICW #*x8000, RO 3 RO/R1 = {XxX 
08 0 SIFD 1 CMPG : compare :X: with 1.0 
iA (18 GEQ_T0_1.0 > branch if !X! >= 1.0 
1 
1 
1 
1 
1 
1 


DOOOCOCOOCOCOS OCOOVGOOO OO OWWOWOWOWWOOOo NNN 


NAME WN 9 OD NOA UE WIN 9 OD NOAUE WIN 8 ODNAU EWN O OONOU EU" 
@o 
7) 
m 
© 


50 50 44FD MULG RO, RO : RO/R1 = Xe2#2 
50 08 50 43FD SUBG RO, #1, RO ; RO/R1 = 1.0 - X¥e2 
og 16 JSB MTHSGSORT_RS + RO/R1 = GSQRT(1-X##2) 
56 O 47FD DIVG3 ~ R R 3; RO/R1 = X/GSQRT(1-X**2) 
00000000 ' GF 17 JMP G*MTHSGATAN_R7 : RO/R1 = GATAN(X/GSQRT(1-X**2)) 
: 1 =< ix! 
GEQ_T0_1.0: 
OE 16 , BGTR ERROR : branch to ERROR if ‘xX! > 1.0 
1 
17 3 
14 ; ix! = 1.0 
! ; 
50 A7 AF 7D 1 MOVQ G_PI_OVER_2, RO 3 RO/R1 = PI/2 
56 53F 3 TSTG Rb 3 test the sign of X 
sr 1 1 GEQ ae : branch if 
50 0 52F0 0 MNEGG RO, RO 3; RO/R1 = =-PI/2 
Z RETURN: RSB : return with result in RO/R1 
4 & : 1< txt, error 
64 5: 
64 6 


A, 
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aT SIN ; G Floating Point Sine routine 16-$ eI 1:24:4 AX/VMS Macro V04-00 2 
O08 ATHSGASIN.R? - Special GASIN routine grSEb= 1386 41 as 13 YRTARTE. eRedaricasiw.mars1 2% Sy 
4 ERROR: PUSH 
6 MOV 


return PC from JSB routine 
ASH 


bf ane INVARGMAT, -(SP) : i € ndition value 

RO : = result = {Peserve operand -0.0 
; gees to s SY ° chanism vector 
; ap! a; H_RO/R $0 error handler 
: +k. y the result. 
; sign wiry and use real user's PC 
3 spendent of CALL vs 


A 

00000000'GF 02 FB | CALLS #2, G*MTHSSSIGNAL 
z ; return - RO restored from CHFSL_MCH_RO/R1 
6 


é RSB 


LY 


T 


rages 
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3 
? é 
? 4 
? & 
? 4 
? 4s 
Ri 
Be 
HS 
a? 
76 § 
76 4 
76 55 
76 2$ 
76 5 
76 «258 
76 59 
76 «= 260 
76 61 
76 $¢ 
76 6 
76 64 
76 65 
76 6 
76 6 
Be 
Bi 
76 71 
76 ? 
76 ? 
76 74 
76 75 
76 7 
76 7? 
76 ? 
1076 7? 
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MTHRTL.SRC 
= Standard G Floating GASIND 


HGASIN.MAR; 1 
-SBTTL MTHSGASIND 


+ 
FUNCTIONAL DESCRIPTION: 


; GASIND = G floating point function 
; GASIND(X) is computed as: 
If X = 0, then GASIND(X) = F 
If X= then GASIND(X) = 90 
If X = ={, then GASIND(X) = -90. 
If 9 < ik; < 1, then GASIND(X) = ATAND(X/SQRT(1-X*#2)). 
If 1 < iXt, error. 


; CALLING SEQUENCE: 


gasind.wg.v = MTHSGASIND(x.rg.r) 
INPUT PARAMETERS: 


LONG = 4 
x = 1 * LONG 


IMPLICIT INPUTS: 
OUTPUT PARAMETERS: 
VALUE: G floating arcsine of the argument 
IMPLICIT OUTPUTS: 
COMPLETION CODES: 
SIDE EFFECTS: 
Signals: MTHS_INVARGMAT if ‘xX! > 1 with reserved operand in RO/R1 
(co90ed to the signal mechanism vector CHFSL_MCH ByRi by LIBSSIGNAL). 


R 
Associated message is: "INVALID ARGUMENT’’. Résulf is reserved 
operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 


: define longword multiplier 
; Contents of x is the argument 


none 


none 
none 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


ENTRY MTHSGASIND, “M<IV, R2, R3, R4, RS, R6, R7> 
3 standard call-by-reference entry 
: disable DV (and FU), enable IV 
MTHSFLAG_JACKET : flag that this is a jacket procedure in 
MOVAB G*MTHSSJACKET_HND, (FP) 
: set handler address to jacket 
3; handler 


SIN 


6 Flo 
ATasGas 


50 (04 BC at 


eine. Pe 


95 
3% 


int Sine routing soup TECREECHSEE 1:34:48} 


AX/VMS Mac 69 v04-00 Page | 
MTHRTL.SRC SATHGASIN.MAR: 1 (6) 


; : § se of an error in gute 


an Stone convert 


; and resign 
0 


= ne = @value(AP) 


; ty spe 
; return w 


i 


ial GASIND routine 
th result in RO/R1 


to user PC 
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SIND_R7 = Special GASIND routine SEP- 9 CMTHRTL.SR THGASIN.MAR; 1 (7) 
4 99 -SBTTL MTHSGASIND_R? = Special GASIND routine 
4 ° : Special GASIND - used by the standard routine and direct JSB call. 
7 : : CALLING SEQUENCE: 
7? 4; save anything needed in RO:R7 
7 53 MOV asee W : input in RO/R1 
7 $ : JSB MTHSGASIND_R7 
i o8 : RSB 3 return with result in RO/R1 
7 9° 
7 10 MTHSGASIND_R7:: : ial GASIND ti 
56 $0 SOD 0087 i HOUG RO, R6 Sowek inna 
1 : 1g BEQL D_RETURN : return DASIND(O) = 0 
D 14; 
a8p 195; 0< ix! 
Ba 
50 8000 8F AA 008D 18 BICW #*x8000, RO : RO/R1 = {Xi 
08 50 5IFD 98 19 CMPG RO, #1 3 compare {Xi with 1.0 
18 096 0 BGEQ 0_GEQ_T0_1.0 : branch if {X! >= 1.0 
098 353 
8 g :0< Ix! < 1.0 
gh 
50 50 44FD 838 6 MULG RO, RO 3 RO/R1 = Xe22 
50 08 50 43FD O009C 7 SUBG RO, #1, RO 3 RO/R1 = 1.0 = X#*2 
90000000" EF 16 O0A1 : JSB MTHSGSORT_RS : RO/R1 = GSQRT(1-x*#2) 
50 6 0 “Te OA7 DIVG3 RO, R R 3: RO/R1 = X/GSQRT(1-X**2) 
00000000 ' GF 1 Boas ? JMP G*ATHSGATAND_R7 : RO/R1 = GATAND(X/GSQRT(1-X**#2)) 
: : 1 =< ix! 
0B $ D_GEQ_T0_1.0: 
BO 14 OBe : BGTR ERROR : branch to ERROR if !X! > 1.0 
0B4 40 
0B4 41; 
4 4g s ix! = 1.0 
Ope Sea 
50 FF58 CF 7D 4 45 MOVQ G_90, RO ; RO/R1 = 90 
56 53F 0B9 46 TSTG Rb ; test the sign of x 
98 1 OBC 4 GEQ D_RETURN : branch if xX > 
50 0 52FD rt MNEGG RO, RO ; RO/R1 = -90 : 
05 ; r4 D_RETURN: RSB 3; return with result in RO/R1 
C 51 
C 52 .END 


L 6 
MTHSGASI : G Floating Point Sine routine 16-SEP-19 1:24:49 YAX/VMS Macro V04-00 P 11 
Symbol. table : S-SEP-1986 11:95:19 EMTMRIL. SREIRTMGASIN.MaR:1 "29° {Jy 
D_GEQ_10_1.0 R 1 
D RET R 1 
ERROR R 1 
GEQ_10_1.0 R 1 
6.96 y R 1 
€rP1_ovER_2 DR : 
LONG- = 00 
MTHSSJACKET_HND ‘teeeeeee xX 4 
MTHS$S1GNAL eecceree 0 
THSGASIN 018 RG = 01 
MTHSGASIND B28 RG 1 
THSGASIND _R7 RG 1 
THSGASIN_R7 B88 9 RG 1 
MTHSGASIN RO 0000029 RG 1 
TANS _R eereeere Xx 0 
MTHSGATAN_R7 eeeeeree =X 00 
MTHSGSQRT_RS eereeeee =X 8600 
MTHSK_INVARGMAT seeeeere =X 00 
R 00000063 R 01 
VALUE = 00000004 
pew ere reer re ewn emacs + 
: Psect synopsis ! 
femacemon amare nae + 
PSECT name Allocation PSECT No. Attributes 
- ABS, 00000000 < 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
-MTHSCODE 000000C3 ¢«¢ 195.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
gm oeeee ee ser eee eee ees eros + 
; Performance indicators ; 
Phase Page faults CPU Time Elapsed Time 
Initialization 31 0:00:00.10 00:00:01.0 
Command processing 112 :00:00.65 9 SS 
pees } sais ; 94 : 3e 8 00:6 t 
ymbol table sor ° :00:00. 
Pass 2 2 ‘% 0:00:08.74 
Symbol table output : . 0:00:00.09 
Psect synopsis output 3 8: : 3: 338 10 
Cross-reference output : ° 00: 09 
Assembler run totals 316 0: . 0:00:12.4 


[he work ing set Limit was 900 pages. 

5179 bytes (11 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 22 non-local and 0 local symbols. 
412 source Lines were read in Pass 1, producing 14 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


MTHSGASIN ; G Floating Point Sine routine = a4 = 1384 91: $4: 43 yAreyes 8 


P 1 
VAX=-11 Macro Run Statistics MTHRTL.S REIN TH Meas TH. MAR; 1 rats 13) 


foes senwen see anes ee ceca scans $ 


! Macro Library statistics ! 


woe tence omen eee are osm omean $ 


Macro Library name Macros defined 
_S255SDUA28: (SYSLIBISTARLET.MLB;2 wow WF 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE =SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGASIN/OBJ=OBJ$:MTHGASIN MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRC 


NT CORPORATION 
D PROPRIETARY 


AH-BT13A-SE 
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